Manifestations of the Roton Mode in Dipolar Bose-Einstein Condensates 



Ryan M. Wilson, Shai Ronen, and John L. Bohn 
JILA and Department of Physics, University of Colorado, Boulder, Colorado 80309-0440, USA 

Han Pu 

Department of Physics and Astronomy, and Rice Quantum Institute, 
Rice University, Houston, Texas 77251-1892, USA 
(Dated: February 8, 2010) 

We investigate the structure of trapped Bose-Einstein condensates (BECs) with long-range 
anisotropic dipolar interactions. We find that a small perturbation in the trapping potential can 
lead to dramatic changes in the condensate's density profile for sufficiently large dipolar interaction 
strengths and trap aspect ratios. By employing perturbation theory, we relate these oscillations to a 
previously-identified "roton-like" mode in dipolar BECs. The same physics is responsible for radial 
density oscillations in vortex states of dipolar BECs that have been predicted previously. 



The study of ultracold atomic and molecular gases is 
notable for its connection to denser condensed matter 
systems. Ultracold gases can show a strong resemblance 
to condensed matter systems such as vortex lattices, su- 
perfluids and Mott-insulators. Part of the attraction 
to these analogies is the ability to control an ultracold 
atomic or molecular gas, allowing researchers to explore 
regions of parameter space that are difficult to access in 
a naturally occurring system. 

A recent example of this connection between ultra- 
cold gases and "conventional" condensed matter systems 
arises in dilute Bose-Einstein condensates (BECs) con- 
sisting of dipolar particles. An early theoretical study 
investigated a gas of dipoles that are free to move in a 
plane, but are confined in the direction orthogonal to the 
plane and that are polarized in this same direction. This 
system is predicted to exhibit an anomalous dispersion 
relation that possesses a minimum at a characteristic mo- 
mentum, reminiscent of the roton dispersion well-known 
in superfluid He Moreover, the depth of this mini- 
mum is controlled by the strength of the dipolar inter- 
action (proportional to the square of the dipole moment 
and the density). If this interaction is large enough, the 
roton minimum can become degenerate with the ground 
state. 

In experiments, however, the gas is confined, leading to 
a discrete excitation spectrum rather than a continuous 
dispersion relation. Nevertheless, signatures of the roton 
excitation have been identified in calculations with fully 
three-dimensional trap geometries Q. Certain excita- 
tions exhibit nodal structures on the same length scale as 
the free rotons. Moreover, the excitation energies of these 
modes drop rapidly as the dipolar interaction strength 
increases. We note that the first dipolar BECs (dBECs) 
have already been created using atomic 52 Cr 3- 6] while 
molecular BECs (promising far larger dipoles and tun- 
able dipole moments) are the target of active experimen- 
tal work. 

In this Letter we explore another aspect of roton 
physics that may be observable in dipolar gases. It has 



previously been suggested that boundaries in superfliud 
4 He, including vortex cores, should give rise to radial 
density oscillations whose length scale is characteristic of 
the roton wavelength 0-0] • More recently, calculations 
of vortex states in a dBEC in a hi ghly oblate trap have 
exhibited similar radial structures [10j, raising the ques- 
tion of the relation between these structures and rotons 
in this system (progress has also been made in the under- 
standing of the vortex state in a dBEC in the Thomas- 
Fermi regime. However, in this regime the vortex does 
not manifest a radial ripple [III EH ) ■ 

Our objective in the present work is to explore this 
relationship. Whereas the complete description of super- 
fluid He is complicated by strong interactions, this is not 
the case in a dBEC, where the gas is dilute enough that 
a mean-field approach should work quite well [13j. In- 
deed, we find that the main effect generating the radial 
density oscillations is the perturbation caused by the cen- 
trifugal potential of the vortex state. This perturbation 
contaminates the ground state with the lowest-lying ex- 
cited state, which, in the limit of strong interactions, is 
the roton excitation. We demonstrate this effect by ap- 
plying a perturbation theory to the nonlinear mean-field 
equations; the perturbative approach is in good agree- 
ment with our full numerical calculations. 

At very low temperatures, N bosons trapped in an 
external potential U(r) may be described within mean 
field theory by the nonlocal Gross-Pitaevskii Equation 
(GPE): 

[_|lv 2 + C/(r) + (iV-l) (1) 
2m 
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where fy(r,t) is the condensate wavefunction (with unit 
norm), r is the distance from the trap center, m is the 
boson's mass, and V(r — r') is the two-particle inter- 
action potential. We consider the case of a cylindri- 
cal harmonic trap, for which the external potential is 
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U(r) = ^muj 2 (p 2 + X 2 z 2 ), where A = uj z /uj p is the t rap 
aspect ratio. The interaction potential has the form [14[ 

A-Kh 2 a s ,9 1 — 3cos 2 , . 

V{T ~ T ) = ~^' 5(r ~ r } + d | r _ r /|3 ( 2 ) 

where a s is the s-wave scattering length, d is the dipole 
moment, and 8 is the angle between the vector r — r' and 
the dipole axis. The first term in V(r — r') is the famil- 
iar contact potential, while the second term is the long- 
range anisotropic dipole-dipole potential. This potential 
describes interactions of dipoles that are polarized along 
the trap axis, as could be achieved in an experiment by 
applying a strong external field. For the sake of illumi- 
nating purely dipolar effects, we set a s = in this work, 



a limit that can potentially be achieved experimentally 
in 52 Cr 0. 

Due to the azimuthal symmetry of both the trapping 
potential and the dipole-dipole potential, the ground 
state solutions of Eq. ([1]) may be written in the form 
^(r, t) = tp(p, z) e lkv , where k is the quantum number 
representing the projection of orbital angular momentum 
about the trap's axis [16j . The k = solutions of Eq. ([T]) 
correspond to rotationless BECs, while the k = 1 solu- 
tions correspond to BECs with singly-quantized vortices. 
The radial structure of the vortex is the same as that of 
a rotationless BEC in a trap with a central potential rep- 
resenting the centrifugal force: indeed, by inserting the 
vortex form written above into Eq. (JT|), one obtains: 
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The centrifugal potential Ti 2 k 2 /2mp 2 is responsible for 
the vortex core (i.e., vanishing density at p = 0). 

Following the systematic mapping of the structure and 
stability of k = dBECs in oblate traps Q , here we un- 
dertake to characterize the structure and stability of a 
k = 1 vortex in dBECs. To characterize the dipolar 
interaction strength we introduce the dimensionless pa- 
rameter D = (N — i where aho = y/h/nujp is the 
radial harmonic oscillator length. The dBEC possesses 
dynamic stability when all of the excited-state Bogoli- 
ubov de Gennes (BdG) eigenenergies are real-valued. As 
was found for k = condensates in reference Q, we find 
that for the vortex state there also exists, at any finite 
aspect ratio, a critical value of D = Z? cr ;t above which 
the k = 1 dBEC is dynamically unstable to small pertur- 
bations, while below it the vortex is dynamically stable. 
We assume that the trap itself is non-rotating, so that 
the vortex is not the lowest energy state and therefore is 
not thermodynamically stable. However, here we are in- 
terested in the question of the dynamical stability, which 
is relevant for a closed system at T = 0. 

We find D cr it for various trap aspect ratios by solv- 
ing the linearized BdG equations, as was done in refer- 
ence [l7| ■ Fig. Q] illustrates the regions of dynamical sta- 
bility for k = 1 and k = dBECs at trap aspect ratios 
that are relevant to this Letter. We concentrate here on 
a specific region in parameter space where we find ripples 
in the vortex density profiles, as illustrated in the insets. 
Whereas previously the ripple was reported for a trap 
with aspect ratio A ~ 100 and required the existence of a 
negative scattering length , we find that vortices with 
ripple structure exist also at milder trap aspect ratios of 
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FIG. 1: The red (thin) dotted line marks the maximum dipole 
strength, for a given trap aspect ratio A, below which a ro- 
tationless (k = 0) dBEC is dynamically stable. The colored 
regions represent the dynamically stable region for a k = 1 
dBEC, while the pink (darker) region is where radial oscil- 
lations with local minima are observed. The inset (a) is an 
iso-density surface plot of a k — dBEC perturbed by a small 
gaussian potential centered on the trap axis, while the inset 
(b) is an iso-density surface plot of a k = 1 dBEC. The pres- 
ence of radial oscillations is clear in both cases. 



A ~ 17 and with purely dipolar interactions (a s = 0). 

It is natural to hypothesize that the appearance of the 
ripple in the vortex structure is related to a roton mode 
which is excited by the centrifugal potential of Eq. ([3]). 
This raises the interesting question, could such a rip- 
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pie also be observed in the ground (non- vortex) state of 
dBEC perturbed by an external potential at the center 
of the trap? Such a perturbation may be realized experi- 
mentally by applying a blue-detuned laser along the trap 
axis, taking the form U'(r) = Aexp (— p 2 /2pfy, where A 
is the height of the Gaussian and po is its width. 

For sufficiently oblate traps, k = dBECs exhibit ra- 
dial density oscillations in the presence of such Gaussian 
potentials. Fig. [2] illustrates the radial profiles of k = 
dBECs in a harmonic trap with aspect ratio A = 17 and 
with a Gaussian potential having A = hoj p and po — 
0.2 Oho- To give a concrete example, for 52 Cr atoms in a 
harmonic trap with radial frequency co p = 2n x 100 Hz, 
this translates to having a beam width of po = 280 nm. 
In this trap, an interaction strength of D — 181.2, very 
near the point on instability for a k — dBEC in a 
trap with the above aspect ratio, may be achieved with 
~ 104,000 52 Cr atoms. It is seen that in this case even 
a small Gaussian perturbation makes a dramatic change 
in the dBEC density profile. The radial oscillations near 
-Dcrit are much more pronounced than for a smaller dBEC 
with D — 100. This is suggestive of the roton presence 
in this structure, since the roton is expected to emerge 
with increasing dipolar interaction strength (he, increas- 
ing density for a fixed dipole moment) [3, [si, [l7[ ■ As was 
shown in reference Q, the roton mode undergoes a sig- 
nificant decrease in energy with increase in D until it 
achieves zero energy at -D C rit, marking the point of dy- 
namical instability for the k — condensate. Beyond this 
_D cr it, the roton energy is purely imaginary. Examining 
the nature of the roton itself within BdG theory tightens 
up its relationship with the observed structure discussed 
above. 

For a k — condensate, the coupled BdG equations 
reduce to a single equation, given by 



GF\f) = w 2 |/). 



(4) 



Here, G = P(G - p)P and F = P(F - p)P, where 
P = I — | V I'){4'| is the projection operator into the 
space orthogonal to ground state wavefuncion Also, 
G = H a + C and F = H + C + 2 X, where H is the 
zero-interaction Hamiltonian, C describes a direct inter- 
action and X describes an exchange interaction between 
the Bogoliubov quasiparticle with eigenvector |/) and the 
condensate. All of these operators are defined as in refer- 
ence [l3|- The eigenvector |/) is given by |/) = \u) + \v), 
where {u, v} are the familiar Bogoliubov eigenfunctions. 
The cj appearing on the right hand side of Eq. (Q} is the 
energy eigenvalue corresponding to |/). 

In Eq. it is understood that the linear space on 
which F and G act, and to which |/) belongs, is orthog- 
onal to \^S). Thus, we eliminate a non-physical solution 



with eigenvalue zero 18(. The justification for working 
in this reduced linear space is that it can be shown that 
all physical excitations obey (f\^) =0 [lij . 
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FIG. 2: Radial profiles of the k — dBEC subject to the 
perturbing potential 17' (r) = huj p exp (— p 2 /2(.2 aho) 2 ) in a 
trap with aspect ratio A = 17. The red dash-dotted line 
represents the trapping potential at z = 0, the black solid 
line represents the radial profile of the dBEC at D — 100 and 
the blue dotted line represents the radial profile at D = 181.2, 
near the point of dynamic instability for the k — dBEC. The 
"+" signs represent the perturbation theory results and the 
thin dotted lines represent the unperturbed radial profiles at 
the corresponding dipole strengths. 



It seems natural to assume that the roton mode dom- 
inates the structure of the perturbed dBEC near insta- 
bility because its energy is much lower than the energies 
of the other BdG modes. To explicitly demonstrate this, 
one needs to formulate a perturbation theory of the non- 
linear GPE with respect to external potential perturba- 
tion. 

To do so, we begin by writing a perturbation to the 
trapping potential as U — > U + U' , where U' is the small 
perturbation. The response of the condensate wavefunc- 
tion to this perturbation is then |^) — > |\E r ) + We in- 
sert these expressions into Eq. (fTJ) , linearize in the primed 
quantities, and obtain the equation 

F|$') = -PU'\$). (5) 

To solve Eq. ©, we introduce a basis defined by the 
eigenvalue equation 
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(6) 



and use its eigenfunction solutions to expand 1^') in the 
\ip n ) basis. Plugging these expansions back into Equation 
(0 and working to first order gives the expression for the 
wavefunction perturbation, 



}^ — : \<Pr. 



(7) 



This derivation involves the use of the orthogonality con- 
dition = and the fact that (</3 n |*) = 0. The 
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FIG. 3: Radial profiles of excitations on a rotationless dBEC 
with dipole strength D — 181.2 in a trap with aspect ratio 
A = 17. The solid blue line represents the BdG roton mode 
while the red marks represent the F-operator eigenfunction 
with eigenvalue fi, \<po). 

final expression is formally identical to that of the usual 
perturbation theory of the linear Schrodinger equation. 

The connection between the BdG roton mode and the 
perturbative modes is clear in the limit that the roton 
mode becomes degenerate with the ground state. In this 
limit, the roton energy oj goes to zero. In Eq. (0}, this 
means that GF has eigenvalue zero. Now, note that the 
operator G is positive semi-definite (its lowest eigenvalue 
is zero, with eigenfunction |^}. This is indeed the ground 
state, since Y$>) is nodeless). Accordingly, the operator 
G that, by definition, acts on the linear space orthogonal 
to \ty), is positive definite. It then follows that any so- 
lution of GF | /roton) = must also satisfy F|/ ro ton) = 0. 
Thus, \ipo) — | /roton) is a solution of Eq. (|6]) with eigen- 
value £o = 0. Since is written as an expansion in 
\ip n ) with weights proportional to l/e n , the eigenfunc- 
tion \<po) with eigenvalue £o ~ makes a contribution 
to 1^') that strongly overwhelms the contributions of the 
other eigenf unctions. Thus, in the limit that the roton 
energy goes to zero, |\J ,/ ) is dominated by the BdG roton 

mode, | /roton)- 

To show that \ipo) becomes identical to BdG roton 
mode |/roton) when the roton energy goes to zero, Fig. [3] 
shows the radial profiles of both of these excited modes 
for a rotationless dBEC with dipole strength D = 181.2 
in a trap with aspect ratio A = 17, which is very near the 
point of instability. Additionally, Fig.[2]illustrates the ac- 
curacy with which this perturbation theory predicts the 
wavefunction of a dBEC when perturbed by a Gaussian 
potential, as discussed earlier in this Letter. 

Recall that the k = 1 solution of the GPE gave rise to 
a centrifugal potential in the radial part of Eq. (jHJl . This 
potential is constant along the trap axis and decreases 
quickly in the radial direction. So, just as the Gaussian 



potential perturbs the dBEC and gives rise to ripples on 
its density profile, we expect similar behavior for trapped 
dBECs with a centrifugal potential, i.e., dBECs with vor- 
tex structure. To treat the centrifugal potential with our 
perturbation theory, we introduce a radial cutoff that is 
chosen to be much smaller than the spatial extent of the 
vortex core itself. We find that for large A there is good 
agreement between our perturbation theory and the re- 
sults of our exact calculations. Just as is the case for a 
Gaussian perturbing potential, the roton mode is respon- 
sible for the rich structure observed in the k = 1 vortex 
state of a dBEC close to instability. 

In conclusion, we have developed a perturbation the- 
ory for the GPE and have applied it to dBECs perturbed 
both by thin gaussian potentials centered on the trap axis 
and centrifugal potentials. This theory allows us to re- 
late the radial oscillations observed on the exact ground 
state profiles of perturbed dBECs to the roton mode ob- 
served in the BdG spectrum of rotationless dBECs. For 
52 Cr and the trap parameters discussed in this Letter, 
the length scale of the oscillations is ~ 2/im. This is in 
comparison to the length scale of the predicted ripple in 
the 4 He vortex, which is of the order of 1 A, and has not 
been resolved experimentally up to now. 
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